Strain localization driven by structural relaxation in sheared amorphous solids 
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A two dimensional amorphous material is modeled as an assembly of mesoscopic elemental pieces 
coupled together to form an elastically coherent structure. Plasticity is introduced as the existence 
of different minima in the energy landscape of the elemental constituents. Upon the application 
of an external strain rate, the material shears through the appearance of elemental slip events 
with quadrupolar symmetry. When the energy landscape of the elemental constituents is kept 
fixed, the slip events distribute uniformly throughout the sample, producing on average a uniform 
deformation. However, when the energy landscape at different spatial positions can be rearranged 
dynamically to account for structural relaxation, the system develops inhomogeneous deformation 
in the form of shear bands at low shear rates, and stick-slip-like motion at the shear bands for the 
lowest shear rates. The origin of strain localization is traced back to a region of negative correlation 
between strain rate and stress, which appears only if structural relaxation is present. The model 
also reproduces other well known effects in the rheology of amorphous materials, as a stress peak in 
a strain rate controlled experiment staring from rest, and the increase of the maximum of this peak 
with sample age. 
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I. INTRODUCTION 



Plastic deformation of crystalline materials is well un- 
derstood in terms of movement of defects (dislocations) 
that are responsible for the energy dissipation that is the 
hallmark of plasticity. On the contrary, the description 
of plasticity of amorphous materials faces the problem 
that there is not a well defined state upon which 'de- 
fects' can be introduced. This has greatly delayed the 
development of a theory of plasticity of amorphous ma- 
terials. For applications, it is clear that such a theory 
would be extremely useful. Many complex liquids [1] and 
polymers Q that do not easily crystallize upon tempera- 
ture reduction, different kinds of foams [3| , and metallic 
glasses Q are prominent examples of amorphous mate- 
rials for which a full theory of yielding and plasticity is 
still lacking. Experiments show that the behavior of such 
materials upon shearing is far from trivial ^] . One of the 
most interesting phenomena occurring when an amor- 
phous system is sheared is strain localization ^ , namely, 
the material deforms plastically only in very narrow re- 
gions called shear bands, where the material is sometimes 
referred to as becoming 'fluid-like' whereas it remains 
solid in the rest. This phenomenon occurs typically at 
sufficiently low shear rates, whereas at higher shear rates 
homogeneous deformation is observed H, 111- Effects 
associated with the time elapsed since the sample was 
prepared (the sample 'age') have also become apparent: 
the static yield stress of these materials is observed to 
increase with age [13, [Hi- Although there are partial 
qualitative justifications to some of these facts, there is 
not a comprehensive description of all the phenomenol- 
ogy 

Probably the most promising theoretical approach to 
the plasticity of amorphous materials at a fundamental 



level is the shear transformation zone (STZ) theory of 
Falk and Langer Within this framework, it is as- 

sumed that elemental pieces of the material can arrange 
in two different structural states. Upon the application 
of an external stress, one of these states is favored with 
respect to the other, and a switching between them may 
occur. These transitions are the elemental plastic events 
in the theory. The STZ theory describes plasticity at an 
almost microscopic level. For instance in simulations us- 
ing bi-disperse systems each STZ is assumed to contain 
two large and two small particles. It seems that the STZ 
theory is able to reproduce some of the phenomenology 
observed in the shearing of amorphous solids, however, 
obtaining full macroscopic predictions seems to be a dif- 
ficult task, requiring the introduction of some ad hoc as- 
sumptions. 

An alternative possibility is to consider plasticity at 
a mesoscopic scale. In this approach (pioneered in (l3j ) 
global plastic effects appear as a consequence of the in- 
terplay between localized plastic events within individual 
constituents of the system, and the elastic distortion they 
produce in the neighborhood. I follow in principle this 
approach, with an important difference in implementa- 
tion (see below). In the absence of other effects, it is 
shown that this leads to an homogeneous (on average) 
deformation upon the application of external strains. An- 
other ingredient is necessary to explain the appearance of 
shear bands. The novel ingredient in the present formal- 
ism is the incorporation of structural relaxation. It will 
be shown through a simplified analysis that this relax- 
ation produces a region of negative dependence of stress 
on strain rate, and then an unstable situation that can 
give rise to shear banding fli^ . These predictions are con- 
firmed with full simulations of the model. In addition, 
I show that structural relaxation is also able to explain 
other phenomena, as a peak in the stress response to an 
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applied strain rate starting from rest. The height of this 

peak is found to increase with sample age. 5C 



II. MODEL 



I model a (two-dimensional) amorphous system at a 
mesoscopic level using the (symmetric) linear strain ten- 
sor eij — {dui/dxj + duj /dxi)/2, with Ui {i — x, y) be- 
ing the components of the two dimensional displacement 
field. For convenience I introduce the three independent 
variables ei = (en + e22)/2, 62 = (en - e22)/2, and 
63 — ^12- In this way, ei is related to volume distortions, 
and 62 and 63 are the two independent shear distortions, 
which are related by a 7r/4 rotation in the x-y plane p^ . 

At the level of description provided by the strain ten- 
sor the system is assumed to be elastically coherent, this 
means that a compatibility condition between the values 
of ei, £2, and £3 must be fulfilled. This is known as the 
Saint Venant condition, that can be written as 

(52 -t- 52)ei - (dl - dl)e2 - 29.5,63 = (1) 

Beyond this global elastic compatibility condition 
(which implies in particular that structural misfits as 
those originated by cracks in the material are not al- 
lowed), the response of the system will be dictated by 
the form of the local free energy /(ei, 62, 63), from which 
the total free energy F is calculated by spatial integra- 
tion. In the case of modeling a perfectly elastic and 
isotropic material we have to choose / to be of the form 



/(ei,e2,e3) 



Be\- 



^^(62-1-63), having a single minimum 



at Ci = and with B and [i related to the bulk and shear 
modulus of the material (see Ref. [l^ for other choices 
appropriate to other physical situations) . To model plas- 
ticity I will take a form for / with minima at a set of 
positions in the plane 62-63. This means that there is 
a set of discrete deformations at which the system is lo- 
cally relaxed. The jump between different minima as the 
system is driven externally by shear will be the elemental 
plastic events of the model. Specifically, / will have the 
form 



/ = /o(e2,e3)+B6?. 



(2) 



The dependence of / on 61 has been assumed to be sim- 
ply quadratic, since plastic effects are associated to shear, 
and not to volume distortions. The function /o describes 
the energy landscape in the 62-63 plane, which should 
have minima in different positions. There are many dif- 
ferent ways of generating such a landscape. The proper- 
ties of the model is not crucially dependent on this choice, 
and the one implemented here is just one possibility. To 
generate a disordered distribution of minima in the 62- 
63 plane (akin to an amorphous system) the following 
scheme is adopted: 



/o(e2,e3) 



A ^ cos \ujQ (62 cos(6') 



■ 63 sin(6')) -I- 09] 
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FIG. 1: A typical form of the local free energy landscape 
given by the function /o, in the plane 62-63. Parameters used 
are: ^Jmin = 0.5, bJraa^ = 1.5, and the sum over 6 in Eq. |3] 
discretized using 13 points. These parameters are maintained 
throughout the paper. 



i.e., /o is chosen as a sum of oscillating terms with differ- 
ent frequencies and phases, in all directions of the 62-63 
plane, parameterized by the angle 9. To keep the numer- 
ical evaluations within reasonable limits the sum over 9 
is discretized using around ten terms. The frequencies 
ujQ and phases (pe of the harmonic components are cho- 
sen randomly (frequencies between ujmin and ujmaxi and 
phases between and 27r) in different spatial positions to 
take into account that the energy landscape at different 
spatial points needs not to be the same. As the locations 
of the minima in the 62-63 plane do not coincide in dif- 
ferent spatial positions of the system, due to the elastic 
compatibility condition ([1]) spatial fluctuations in the lo- 
cal stresses will typically be present. An example of a 
typical form of the function /o is shown in Fig. [T] [T7| . 
Defining the local principal stresses Ui as 



(^i[x,y) 



SF 



5ei{x,yy 



(4) 



the dynamical evolution of the strain is obtain through 
a first order temporal evolution equation of the form 



dei{x,y) 
dt 



Wt{x,y) + Ai{x,y,ei,t) 



(5) 



(3) 



where A^ is a Lagrange multiplier chosen to enforce the 
compatibility condition ^ [15], and 77 sets the time scale. 
In equilibrium {dei(x,y) / dt — 0), this equation reduces 
to the standard elastic equilibrium equations, namely 
d/dxi (SF/Seij) = [l^. Up to here, except for some 
differences in implementation, the model may be consid- 
ered to be rather similar to some previous approaches 

MM. 
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Structural relaxation, which is a key ingredient in the 
model, is introduced as an additional mechanism of free 
energy minimization in the system. In the present model, 
structural relaxation is a relative shift of the local energy 
landscapes at different spatial points, that tends to re- 
duce the total free energy. 

To obtain this effect, a couple of additional fields 
62(0;, y), and e'^{x,y) are introduced in order to get a 
rigid relative displacement of the energy landscapes in 
different spatial positions, by changing the function /q 
according to 



/o(e2,e3) -> /o(e2 



62,63 



4) 



(6) 



The tendency toward more relaxed configurations will be 
modeled by the dynamical evolution of 62 and 63. In the 
simplest possible approach this evolution will be chosen 
to be 



dt 



SF 



de°{x,y) 



(7) 



(« = 2,3). Note that a conserving dynamics (i.e., the 
inclusion of the Laplacian operator) is used to guaran- 
tee that uniform stresses do not relax. The value of the 
coefficient A measures the intensity of structural relax- 
ation. Although the dynamics of e^'s and e"'s look al- 
most identical (beyond the fact the one is non-conserving 
and the other conserving) there is an important differ- 
ence between and e° in the fact that the e^'s satisfy a 
compatibility constrain, whereas e°'s do not. Also, the 
dynamics of the (that models structural relaxation) 
will be assumed to be much slower than that of the e^'s, 
that is governed by elasticity. In a situation in which ex- 
ternally applied strains are fixed in time and uniform in 
space, evolution through ^ and ([7]) produces eventually 
a configuration in which the principal stresses ai (and 
also the usual stress tensor of the system) become uni- 
form throughout the system. In this state the system has 
reached its most relaxed configuration locally available. 



III. RESULTS 

To gain insight into the kind of response we will observe 
in the system, it is good to start with a simplified analysis 
of the model equations. This is done in the next sub- 
section, and the results with the full model are presented 
after that. 



A. Simplified Mean Field Analysis 

Let us consider an extremely simplified version of the 
model, that will be useful to understand roughly the kind 
of results we will obtain in the full version. I want to em- 
phasize that this 'simplified version' is not meant to be 



a controlled approximation on the main model in a par- 
ticular parameter region, but just a simpler, 'similarly 
looking' model, that admits a more transparent analysis 
of the effect of structural relaxation. The simplified ver- 
sion considers only one of the shear strains, the one that 
is driven externally by the applied strain. Let us call this 
strain generically e (in particular cases this could be 62 or 
63). It is assumed that the other components of strains 
adjust to satisfy the compatibility condition, in such a 
way that we do not need to take it into account explic- 
itly. The external control parameter will be the mean 
strain in the system e. I assume e changes linearly in 
time, i.e., e = jt, 7 defining the strain rate. Aiming at 
a mean field description of this simplified model, I will 
consider a dynamical evolution for e{x,y) of the form 



de{x,y) 
dt 



S f{e{x,y)) 
Se{x,y) 



+ hie-eix,y)) (8) 



The last term takes into account the coupling between 
the values of e in different points of the system through 
a mean field interaction. The local free energy / has 
minima at different values of e. To have a simple form 
for / I take 



f{e{x,y)) = Asm[uj{x,y) (e - (p{x,y))] 



(9) 



i.e., a simple harmonic form, with amplitude and phase 
that vary across the system. Structural relaxation is in- 
cluded by the dynamical evolution of the phases (f>{x,y), 
again in a mean field way, through 



d(t){x,y) 
dt 



= -A 



5f{e{x,y)) S~f 
5(p{x,y) 5^ 



(10) 



Considering first the case of no structural relaxation 
(A = 0), it can be noticed that the present model reduces 
to that considered by Fisher [l^ to study sliding charge 
density waves, in mean field approximation. Using the 
results in Ref. [l^, we can conclude that in the absence 
of structural relaxation and in the strong pinning regime 
[h/{Arii) sufficiently small, I assume this is the case], the 
stress-strain rate curve for the model starts at a finite 
value of stress, and from this value it increases with ap- 
plied strain rate 7 as 7^^'^, i.e, the stress-strain rate curve 
satisfies, for low strain rates an equation of the form 



(7(7) (To + Cj 



2/3 



(11) 



The main interest in analyzing the present model is to 
consider the effect of finite structural relaxation. We see, 
according to (fTU]) , that relaxation processes tend to alter 
the phases 4> throughout the system and reduce further 
the free energy. If we think of a situation in which all the 
phases have become equal, we understand clearly that 
the system will be able to support larger static stresses 
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FIG. 2: Points: Stress vs strain rate results for the simplified 
model in mean field approximation {A=l, r] = 1, h = 0.5, and 
the frequencies in Eq. ([9]) chosen randomly in the interval 0.5- 
1.5). Lines: Fitting to a function of the form given by Eq. 
(|12[) . with a single set of parameters ctq, C and D. Values 
of A are indicated. For no structural relaxation A = 0, the 
numerical results follow closely the analytical result obtained 
by Fisher ([!§], and Eq. JTlJ). In the presence of relaxation 
A 0, an additional contribution of the form suggested by 
Eq. [T2] is obtained. Note in particular that now the curve 
gets a regime of increasing stress upon strain rate reduction. 

that in the absence of relaxation, because now different 
parts of the system support the applied strain coherently. 
On the other hand, the temporal scale of shearing and re- 
laxational process (whose ratio is controlled by the quan- 
tity 77/ A) will indicate whether relaxational effects have 
enough time to influence the dynamics. An order of mag- 
nitude estimation of the effect of relaxational terms indi- 
cates that the u vs 7 curve for the whole system will get 
a correction more and more important as 7 is lower and 
lower. This estimation leads to 

a(7) ^ao + C^'"^ + BA/(r,7) (12) 

This estimation should be considered only as a qual- 
itative rule of thumb. In particular, the divergence for 
7 — > is spurious. But it already manifests the crucial ef- 
fect of relaxation: a can get a regime in which it increases 
as 7 is reduced. To confirm this estimation, in Fig. [5] I 
present numerical results for this simple model, with and 
without structural relaxation. For 7 = the analytical 
result of Fisher is reproduced. For 7 ^ we get in fact a 
regime in which a increases upon 7 reduction. The trend 
is well described by a I/7 behavior. 

The behavior we have found in this mean field treat- 
ment of this toy model will produce important conse- 
quences in the more complete simulation presented in 
the next section. In fact, a non-monotonic form of this 
curve has been suggested to be crucial in the appear- 
ance of shear bands in some phenomenological models 
[1^ . In a certain sense, we will be able to consider 
the model in the next section as an assembly of pieces 
with a non-monotonous stress-strain rate behavior, as the 



one obtained here. Two salient features can be already 
described. It is easy to see that when stress increases 
upon strain rate reduction, a uniform deformation rate 
throughout the system will be an unstable situation: if 
some region of the system shears more rapidly than other 
region, stress diminishes in the former and increases in 
the latter, leading to the result that shear is completely 
halted in some part of the system and becomes spatially 
localized. The elastic compatibility in the system pro- 
duces that the regions where shear is accumulated take 
the form of bands. This is the origin of shear bands, that 
will be observed in the simulation of the next section. A 
second consequence of non-monotonic stress-strain rate 
behavior will be the appearance of stick-slip like motion, 
as it is well known for instance in models of sliding be- 
havior, where typically the 'static' friction coefficient is 
larger than the 'dynamical' one (20j . 

B. Full Numerical Solution 

On the basis of the results presented in the previous 
section, in particular, on the effect of structural relax- 
ation, I will present now the results of full numerical 
simulation of the model. The results presented were 
obtained in systems of size 128x128 lattice points with 
periodic boundary conditions, upon the constrains that 
ei = 63 = 0, 62 = e, where bar notes spatial averages 
over the system. Then e is the main control parame- 
ter, which represents a shear strain (I will call it simply 
'strain', from now on) that compresses the system in the 
y direction and expands it in the x direction, leaving the 
volume unchanged. Note that this distortion will show 
a tendency to produce shear bands at 45 deg of the x-y 
axis. The reported stress g is defined as cr = 02- Di- 
mensionless units are used for all quantities by choosing 
5 = 1.5,77 = 1,^ = 1. When necessary, the share rate 
7 is defined, in such a way that e = ^t. 



1. Results for uniform strain rate conditions 

As pointed out in the previous section, the stability or 
instability of a homogeneous deformation in a system will 
be related to a monotonic or non-monotonic form of the 
stress vs. strain rate curve. Thus as a starting point, in 
Fig. [3] I present the results of this calculation for differ- 
ent degrees of structural relaxation. In the curve without 
relaxation (A = 0) we see a monotonous increasing of a 
with 7. A rather good fitting is obtained by an expres- 
sion of the form a — gq + Cj^^'^ . A closer examination 
of the deformation mechanisms in real space, shows that 
deformation proceeds through a sequence of elementary 
plastic events. They are generated by the abrupt jump 
of an elemental piece of the system from one minimum to 
another one. This change generates a quadrupolar elas- 
tic field that disturbs the neighborhood, and can trigger 
additional transitions in nearby points. These slip events 
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FIG. 3: Stress vs strain rate curve of the complete model. 
The results show the same trends seen in the previous figure 
for the simpler model. 



do not span in general the whole sample, but are local- 
ized in space. In Fig. [H^a) we see a snapshot of the 
accumulated deformation during a rather small period of 
time, in which the existence of some of these events is 
apparent. However, when longer periods of time are con- 
sidered, it is found that the elementary slip events occur 
through the whole system [Fig. IDJc)], and on average, 
the deformation becomes homogeneous. 

To obtain shear bands, i.e., deformations that remain 
localized for long periods of time, structural relaxation 
must be included. First of all, we see in Fig. [3] that now 
the stress vs. strain rate curve becomes non-monotonic, 
as it was the case in the simplified analysis of the previ- 
ous section [2ll |. The deformation mechanism in this case, 
corresponds also to a sequence of elemental slip events, as 
in the previous case, but now, when deformation is accu- 
mulated over long periods of time, its spatial distribution 
shows the appearance of a shear band, where deforma- 
tion is localized. This can be seen in Fig. [S] where I 
show the accumulated deformation for different values of 
strain rate. While deformation remains uniform at rel- 
atively large strain rate, it develops for sufficiently low 
strain rate a clear localized form consisting of a shear 
band in which the system is sheared, and the rest of the 
system remains in a static configuration. 

The importance of structural relaxation in the appear- 
ance of shear bands becomes clearer by plotting the spa- 
tial distribution of stresses in the system (Fig. [6]). This 
is precisely the variable that tends to be homogenized by 
the structural relaxation terms. At large shear rates [Fig. 
[SJa)], when deformation is uniform, we see that a keeps 
large fluctuations in the whole system. This is due to 
the large value of shear rate, that prevents the structural 
relaxation to act efhciently. However, at lower shear rate 
we see that the region that is not shearing has acquired 
a rather uniform distribution of stresses whereas the re- 
gion that is shearing keeps larger stress fluctuations. This 
fact makes the non-shearing region even stronger against 
shear, which is then forced to remain localized in the 




FIG. 4: Deviation from average deformation in a simulation 
without structural relaxation, for progressively larger peri- 
ods of time to ((a) to = 10, (a) to = 10^, (a) to = 10^), 
with a fixed strain rate 7 = 0.005. The short time struc- 
ture shows the existence of discrete slip events. They have a 
typical quadrupolar symmetry, and can trigger similar events 
along the diagonals. Over long periods of time, all the sam- 
ple participates of the deformation, that becomes spatially 
uniform. 



shear band. 

The appearance of a shear band can be considered as a 
phenomenon analog to phase coexistence at a first order 
phase transition. This situation occurs here because of 
the non-monotonic behavior of stress vs strain rate curve. 
In fact, in a strain rate controlled experiment, at large 
values of this parameter, the deformation distribution is 
uniform. As 7 is reduced, we reach the reentrant region 
of the curve (the coexistence region), and now part or 
the system adopts a configuration with zero deformation 



6 




FIG. 5: Accumulated deformation in simula- 

tions at progressively lower strain rates, namely 
7 = 0.005,0.001,0.0005,0.0001 from (a) to (d), in the 
presence of structural relaxation (A = 0.001). The total 
average deformation is the same in all cases. Gray scale 
goes from white for no deformation to black for maximum 
deformation. The direction z along which the profiles in Fig. 
[7] are calculated is indicated. 



rate, whereas the rest retains its original velocity. Upon 
a further reduction of 7, a larger fraction of the system 
goes to a state of zero strain rate, i.e, the shear band 
becomes progressively thinner. This argument suggests 
that in the shear banding regime, the deformation rate 
at the shear band is independent of the global deforma- 
tion rate, it is the width of the shear band what changes 
with 7. As a verification of this behavior, in Fig. [7] we 
see profiles of the local deformation rate across the di- 
rection along which shear band forms. We see in fact 
that when the shear band forms, it tends to maintain 
within itself a rather constant value of deformation rate. 
If the global deformation rate is reduced, the shear band 
becomes thinner. There is however a limit to this shear 
band thinning process upon deformation rate reduction 
that allows to understand the origin of another very inter- 
esting phenomenon. The width of a shear band cannot 
be lower than some value, which in the present model 
is related to the size of simulation cells, and in experi- 
mental situations has been observed to be, in granular 
systems, of the order of ten particle diameters [l^]. If 
the deformation rate is too low, even the thinnest shear 
band has to reduce its deformation rate to satisfy the ex- 
ternal constrain. Since a lower deformation rate at the 
shear band would situate it in the negative slope part 
of the stress-strain rate curve, this situation would be 
unstable. The system adopts a dynamics that resembles 
very much the 'stick-slip' motion occurring for instance 
in many examples of friction motion (20| . The system 
alternates between periods in which it rigidly follows the 
external deformation, and abrupt bursts at which slip 




FIG. 6: Instantaneous spatial distribution of the stress a in 
the system, corresponding to the snapshots in the previous fig- 
ure. We see that due to structural relaxation, stress becomes 
uniform in regions that remain rigid, while large fiuctuations 
remain within the shear bands. 

at the shear band occurs. The fluctuations in the in- 
stantaneous stress, and its temporal evolution is illus- 
trated in Fig. [51 where we see in fact the appearance 
of large fluctuations when we enter the regime of a very 
thin shear band. This transition however does not seem 
to be abrupt but smooth. 



2,. Effects of sample preparation 

In addition to the steady state properties discussed in 
the previous section, other interesting properties of amor- 
phous materials have to do with the response to different 
preparation conditions. They can be also discussed in 
the present model. 

Fig. [HI shows the response of the system in simula- 
tions done with two different protocols. In Fig. [H^a) I 
show results for the evolution of stress as a function of 
strain, for a fixed value of 7. The system is started in 
a configuration with ei{x,y) = 0, and it is allowed to 
relax (using a finite A) during some 'waiting time' be- 
fore the strain rate is applied. After this initialization, 
structural relaxation is turned off (A = 0) and a finite 
7 is applied. For sufficiently small initial relaxation the 
curves grow monotonously toward a saturation value at 
large deformation. However, in the presence of an initial 
relaxation, a peak in the transient response is observed. 
This is understood as originated in the more relaxed orig- 
inal structure that is being simulated in this case. Once 
the system has escaped from the deeper free energy min- 
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FIG. 7: (a) Average strain rate across the shear band, for pro- 
gressively lower values of 7, in simulations with a relaxation 
parameter A = 5 x 10~*. The total length plotted corresponds 
to half of the diagonal in Fig. (5) (b) maximum strain rate in 
the system, as a function of the average strain rate. Each dot 
corresponds to one curve in (a). 
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FIG. 8: Standard deviation ha of the instantaneous stress in 
the system, for A = (lower curve) and A = 5 x 10~* (upper 
curve), and examples of the actual evolution at particular 
points (each inset panel depicts the evolution during a time 
interval of = 5). A non-zero A produces an appreciable 
increase in fluctuation at low strain rates, and the evolution 
of a acquires a stick-slip like behavior. 



FIG. 9: (a) Stress vs. strain curve for samples allowed to 
age during the indicated waiting time in the presence of 
structural relaxation (A = 0.01). Shear rate is 7 = 0.01. 
Structural relaxation was turned off when strain is increased. 
Inset shows the the evolution of the height of the peak with 
calculated by first subtracting the results for — 0. (b) 
Response of an aged sample (t^ — > 00) for different shear 
rates, as indicated. Inset shows the height of the peak with 
respect to the asymptotic value. (For this figure a system of 
size 256 x 256 was used) 



imum in which it was located, the disordered nature of 
the minima in different spatial points in the system re- 
stores a totally disordered distribution of minima, and 
stress stabilizes at a value independent of the original re- 
laxation. We can say that after many plastic events the 
system looses memory of how it was prepared. 

In Fig. m^b), we see results for a system allowed to 
relax completely before strain rate is applied (t^, —^ 00). 
Different curves correspond to different values of strain 
rate. As 7 increases, the height of the peak and the 
asymptotic shear stress both increase. This behaviors 
agrees qualitatively with simulations in binary Lennard- 
Jones fluids §1. 



IV. CONCLUSIONS 

In summary, I have adapted an elastic model previ- 
ously used in other contexts [isj to study the appearance 
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of strain localization during the shearing of amorphous 
materials. The main ingredients of the model are: long 
range elastic interactions that are enforced through an 
elastic compatibility condition; a set of minima in the 
shear plane 62-63 available for each discrete element in 
the system; and a mechanism for structural relaxation 
that tends to make the stresses uniform through the sys- 
tem. It is in this last ingredient where the novelty of the 
present approach resides. I have obtained an important 
number of results that have been observed in the shearing 
of amorphous materials, namely, the increase of stress as 



a function of strain rate, the appearance of shear bands 
at low shear rates, and the stick-slip-like motion at very 
small shear rates. Aging effects in the value of the peak 
stress stress of the model are also obtained. The sim- 
plicity of the present approach, which is based on quite 
natural assumptions makes it a useful tool to investigate 
the rheological properties of amorphous materials. 
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